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Complex networks have acquired a great popularity in recent years, since the graph representation 
of many natural, social and technological systems is often very helpful to characterize and model 
their phenomenology. Additionally, the mathematical tools of statistical physics have proven to be 
particularly suitable for studying and understanding complex networks. Nevertheless, an important 
obstacle to this theoretical approach is still represented by the difficulties to draw parallelisms 
between network science and more traditional aspects of statistical physics. In this paper, we explore 
the relation between complex networks and a well known topic of statistical physics: renormalization. 
A general method to analyze renormalization flows of complex networks is introduced. The method 
can be applied to study any suitable renormalization transformation. Finite-size scaling can be 
performed on computer-generated networks in order to classify them in universality classes. We also 
present applications of the method on real networks. 
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I. INTRODUCTION 

Many real systems in nature, society and technology 
can be represented as complex networks [ll, H, H, 0, B S] ■ 
Independently of their natural, social or technological ori- 
gin, most networks share common topological features, 
like the "small-world" property Q and a strong topolog- 
ical heterogeneity. The small-world property expresses 
the fact that the average distance between nodes, as de- 
fined in the graph-theoretical sense, is small with respect 
to the number of nodes, and typically grows only loga- 
rithmically with it. Networks are topologically heteroge- 
neous in that the distributions of the number of neigh- 
bors (degree) of a node are broad, typically spanning 
several orders of magnitude, with tails that can often be 
described by power laws (hence the name "scale-free net- 
works" [8]). 

While "scale-free-ness" implies the absence of a char- 
acteristic scale for the degree of a node, it is not a 
priori clear how this can be related to the notion of 
self-similarity, often studied in statistical physics, and 
also typically related to the occurrence of power laws. 
In this context, several recent works have focused on 
defining and studying the concept of self-similarity for 
networks. The notion of self-similarity is related to 
a renormalization transformation, properly adapted to 
graphs, introduced by Song et al. Q. The renormal- 
ization procedure is analogous to standard length-scale 
transformations, used in classical systems [lol lUl, and 
can be simply performed by using a box covering tech- 
nique interpreted in a graph-theoretical sense. The anal- 
ysis of this transformation in real networks [§] has re- 
vealed that some of them, such as the WWW, social. 
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metabolic and protein-protein interaction networks, ap- 
pear to be self-similar while others, like the Internet, do 
not. Self-similarity here means that the statistical fea- 
tures of a network remain unchanged after the applica- 
tion of the renormalization transformation. Many succes- 
sive papers have focused on this subject, performing the 
same analysis on several networks, introducing new box- 
covering techniques and trying to explain the topological 
differences between self-similar and non-self-similar net- 
works [H, |il,M E K 113, IHl (for a recent review on 
this topic see |19|). 

In this context, the analysis of renormalization flows 
of complex networks |20i] represents a new perspective 
to study block transformations in graphs. Differently 
from all former studies, the study of Ref. [20] is not de- 
voted to observe the effect of a single transformation, 
but to analyze the renormalization flow produced by re- 
peated iterations of the transformation. Starting from 
an initial network, the iteration of the renormalization 
procedure allows to explore the space of network con- 
figurations just as standard renormalization is used to 
explore the phase space of classical systems in statistical 
physics [i^jllH- For these reasons, the analysis of renor- 
malization flows of complex networks represents not only 
an important theoretical step towards the understanding 
of block-transformations in graphs, but also a further at- 
tempt to link traditional statistical physics and network 
science. 

In this paper, we substantially extend the analysis pre- 
sented in [20| . We perform a numerical study of renor- 
malization flows for several computer-generated and real 
networks. The numerical method is applied to different 
renormalization transformations. For a particular class 
of transformations, we find that the renormalization flow 
leads non-self-similar networks to a condensation tran- 
sition, where a few nodes attract a large fraction of all 
links. The main result of the paper lies in the robustness 
of the scaling rules governing the renormalization flow 
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of a network: independently of the transformation, the 
renormalization flow of non-self-similar networks is char- 
acterized by the same set of scaling exponents, which 
identify a unique universality class. In contrast, the 
renormalization flow of self-similar networks allows to 
classify these networks in different universality classes, 
characterized by a set of different scaling exponents. 

The paper is organized in the following way. In sec- 
tion |TT1 we describe the standard technique used in order 
to renormalize a network and define the renormalization 
flow of a graph. We then start with the analysis of renor- 
malization flows of different networks. In the case of 
computer-generated graphs, we distinguish the behavior 
of non-self-similar (section IIII A|) and self-similar (sec- 
tion |IIIB]) networks. Section Hvl is devoted to the analy- 
sis of the renormalization flows of real complex networks. 
Finally, in section |V] we summarize and comment the re- 
sults. 



II. RENORMALIZING COMPLEX NETWORKS 

Differently from classical systems, graphs are not em- 
bedded in Euclidean space. As a consequence, standard 
length-scale transformations cannot be performed on net- 
works since measures of length have a meaning only in 
a graph-theoretical sense: the length of a path is given 
by the number of edges which compose the path; the dis- 
tance between two nodes is given by the length of the (or 
one of the) shortest path(s) connecting the two nodes. 
Based on this metrics. Song et al. 0] proposed an orig- 
inal technique for renormalizing networks (see Fig. [T]). 
Given the length of the transformation , their method 
is given by the following steps: 

• Tile the network with the minimal number of boxes 
Nb; each box should satisfy the condition that all 
pairs of nodes within the box have distance less 
than £b- 

• Replace each box with all nodes and mutual edges 
inside with a supernode. 

• Construct the renormalized network composed of 
all supernodes: two supernodes are connected if 
in the original network there is at least one link 
connecting nodes of their corresponding boxes. 

The former recipe represents a transformation ap- 
plicable to any unweighted and undirected network lead- 
ing to the generation of a new unweighted and undirected 
network, the "renormalized" version of the original one. 
In principle, there are many ways to tile a network and 
therefore the transformation R^^ is not invertible. More- 
over, finding the best coverage of a network (i.e., the one 
which minimizes the number of boxes Nb) is computa- 
tionally expensive: up to now, the best algorithm intro- 
duced in this context is the greedy coloring algorithm [l^ 
(GCA), a greedy technique inspired by the mapping of 
the problem of tiling a network to node-coloring, a well 




FIG. 1: Renormalization procedure applied to a simple graph, 
(left) The original network is divided into boxes and the renor- 
malized graph (right) is generated according to this tiling. 

known problem in graph theory plj . An analogous tech- 
nique, leading to a qualitatively and quantitatively sim- 
ilar transformation Rrg, is random burning (RB) p^ . 
In RB boxes are spheres of radius tb centered at some 
seed nodes, so that the maximal distance between any 
two nodes within a box does not exceed 2rB- Nodes in 
boxes defined through the transformation i^^B satisfy the 
condition defining boxes of the transformation Rg^ , for 
£b = '2.rB + 1. However, the search for minimal box cov- 
erage is much more effective for the GCA than for RB, 
and this may occasionally yield different results, as we 
shall see. 

The strict meaning of self-similarity is that any part of 
an object, however small, looks like the whole [2^]. Sim- 
ilarly, complex networks are self-similar if their statisti- 
cal properties are invariant under a proper renormaliza- 
tion transformation. Song et al. Q have shown that the 
degree distribution of several real networks remains un- 
changed if a few iterations of the renormalization trans- 
formation are performed. Moreover, when this feature 
is verified, the number of boxes Nb needed to tile the 
network for a given value of the length parameter Ib de- 
creases as a power law as £b increases: 

NB{iB)-iB'^ . (1) 

The exponent ds is called, in analogy with classical sys- 
tems, the fractal exponent of the network This prop- 
erty has been verified for several real networks in various 
studies [1, [13, On the other hand, not all real net- 
works are self-similar, i.e. Eq.((T]) and the invariance of 
the degree distribution do not hold for them. For consis- 
tency, these networks are called non-self-similar. 

In contrast to former studies which mostly dealt with 
a single step of renormalization, we are interested here in 
analyzing renormalization flows of complex networks, i.e. 
the outcome of repeated iterations of the renormalization 
procedure described above. Starting from a graph Go, 
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with A^o nodes and Eq edges, we indicate as Gt (with 
Nt nodes and Et edges) the network obtained after t 
iterations of the transformation R: 

Gi ~ R (Go) , G2 — R (Gi) = R^ (Go) , . . . 

... , Gt^RiGt-i) = ... = R'iGo) ■ ^ ' 

Note that in Eq.([2]) we have suppressed the subscript £3 
(or rs) for clarity of notation. In our analysis, we follow 
the flow by considering several observables. We mainly 
focus on the variables 



graphs fSec. lIIIB| . On the other hand, we obtain x* = 
in all cases, except for the particular transformations ob- 
tained for = 1 and £b = 2 on non-self-similar networks 
fSec. IIII A ip . In the next sections, we show our numeri- 
cal results, obtained from the analysis of renormalization 
flows of computer-generated networks, distinguishing be- 
tween the various cases. All values of v and x* are listed 
in Table [J We emphasize the importance of the fact 
that the exponent v is able to classify artiflcial networks 
in different universality classes. 



K,^Kt/{Nt^l) , (3) 
where Kt is the largest degree of the graph Gt , and 



(4) 



which is basically the average degree of the graph Gt 
divided by two. Kt and ijt can assume non-trivial values 
in any graph, excluding trees (for which r]t — 1, Vi). We 
monitor also the fluctuations of the variable Kt along the 
flow by measuring the susceptibility 



Xt = A^o ((^?) - i^t)') 



(5) 



where (•) denotes averages taken over different realiza- 
tions of the covering algorithm. Moreover, we consider 
other quantities like the average clustering coefRcient 
Gt All these observables are monitored as a func- 
tion of the relative network size xt = Nt/No, which is a 
natural way of following the renormalization flow of the 
variables under study. 



III. COMPUTER-GENERATED NETWORKS 

We first consider artificial networks. In the case of 
computer-generated networks, it is in fact possible to con- 
trol the size Nq of the initial graph Gq and to perform the 
well known finite-size scaling analysis for the renormal- 
ization flow. For every computer-generated graph and 
every transformation R, we find that the observable Kt 
obeys a relation of the type 



Kt^F 



{xt-x*)N^ 



(6) 



where -F (■) is a suitable function depending on the 
starting network and the particular transformation used. 
Analogous scaling relations hold for the other observ- 
ables {j]t and Gt) we considered. The susceptibility Xt 
needs an additional exponent 7 since it obeys a relation 



of the type: Xt = K^'G 



with G(-) 



{xt-x*)N'^' 

suitable scaling function. In general, the scaling expo- 
nent V does not depend on the particular transformation 
R used to renormalize the network, but depends on the 
starting network Go: we always obtain v = 2 for any 
non-self-similar network (Sec. IIII A[) and values of v de- 
pending on the initial network in the case of self-similar 



A. Non-self-similar networks 

We consider several computer-generated networks for 
which Eq.([T|) does not hold. Eq.® is able to describe 
the renormalization flows of any of these network models. 
The scaling exponent v = 2 identifies a single universal- 
ity class for all these models. The only difference is given 
by the finite value of > obtained when renormaliza- 
tion is performed with = 2 or = 1 (Sec. IIII A II) . 
Instead, for £b > 2 and rs > 1 we always obtain x* = 
fSec. lIIIA2|) . 



1. rs = 1, £3 = 2. 

For rs — I or £b — the transformation R has a par- 
ticular behavior. In the case of GCA and £3 = 2, at each 
stage of the renormalization flow, the boxes in which the 
network is tiled have the peculiarity to be fully connected 
subgraphs or cliques [231] . In the case of renormalization 
performed with RB and rs = I, spheres are compact sub- 
graphs composed only of neighbors of the selected seed 
nodes. In both cases, at each stage of the renormaliza- 
tion flow, the contraction of the network is much slower if 
compared with the same transformations run for higher 
values of or . 

In Fig. [21 we show some numerical results obtained 
following the renormalization flow of the Erdos-Renyi 
(ER) 24] model with average degree (k) = 2. For both al- 
gorithms used for renormalizing the networks, we clearly 
see a point of intersection between all the curves occur- 
ring at a particular value x* > 0. Interestingly, the same 
values of v and x* hold also for the susceptibility Xt and 
the average clustering coefficient Gt. Numerical results 
for both quantities and their relative scaling are reported 
in Fig. [3 

The same behavior is observed for all the non-self- 
similar networks that we have studied. To mention a 
few, we performed numerical simulations also on the 
Barabasi-Albert (BA) model 8] and its generalization 
given by scale-free networks generated with linear prefer- 
ential attachment [2^. We report in Fig.[4lthe numerical 
results obtained for the BA model: the quantities Kt and 
rjt are shown as a function of the renormalization flow's 
variable xt. Again, a clear crossing point a;* > can be 
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FIG. 2: Study of renormalizatioii flows on the ER model with (k) = 2. The box covering has been performed by using RB with 
rs = 1 (a) and GCA with £b = 2 (b). The figures display m (a,b top) and rjt (a,b bottom) as a function of the relative network 
size Xt- The insets display the scaling function of the variable {xt — x*)Ng^'^ for Kt and -qt- Here the scaling exponent u — 2 in 
both cases. Note that the flow of the renormalization procedure goes from larger (right on the x-axis) to smaller values (left) 
of Xt. 



seen in this case. More importantly, both variables Kt 
and r]t obey Eq.® with = 2. 

The existence of a non-vanishing x* is a peculiarity of 
the renormalization obtained for £b = 2 and rs = 1: 
X* > implies the existence of a special stable fixed 
point, which holds in the limit of infinite network size. 
The fixed point is reached in a number of iterations which 
scales logarithmically with the initial size of the network, 
while the number of renormalization stages needed to 
reach any xt < x* diverges almost linearly with the ini- 
tial system size (see Fig.[S]). Interestingly, the fixed point 
statistically corresponds to the same topological struc- 
ture, independently of the topology of the initial net- 
work (see Fig. ^ . This particular fixed point is a graph 
where a few nodes attract a large fraction of all links 
[i.e., Kt (x*) w 1]; such hub nodes have degrees which are 
distributed according to a power law (see Fig. [6^) . More- 
over, the network obtained at the fixed point is composed 
of nodes with clustering coefficient (C) and average de- 
gree of the neighbors (fc„„) which decrease as a power law 
as the degree of the node increases (see Figs.[6)D and[6t, 
respectively). Fig.[5]i is a graphical representation of the 
graph obtained at x* when starting from an ER network 
with A'o = 30000. The presence of many star-like struc- 
tures gives an explanation of the results described above 
(Fig. El): for — 1, the center of the first chosen box will 
be with high probability a low-degree node ( "leaf" in the 
figure), so the box will include only a low-degree node 
and the attached hub, and the other low-degree nodes 
(the other leaves of the star) will need one box each. This 
is why Nt decreases very slowly. Renormalization steps 
with 'small' boxes (r^ = 1 or = 2) make it therefore 



'difficult' or 'slow' to modify appreciably such structures. 



2. rs >l, Ib> 2. 

Renormalization fiows with £b > 2 and > 1 have 
been discussed in [20], and we present some additional 
results in Fig. [T] The renormalization flows of non-self- 
similar computer-generated networks still obey Eq.® 
with V = 2, and the main difference with respect to the 
particular cases £b = 2 and tb = ^ consists in the value 
of X* , which is now x* = 0. As usual in statistical physics, 
however, the precise value of the threshold is less relevant 
than the value of the exponents describing the flow. In 
this respect, the robustness of the value v = 2 strikingly 
shows that all non-self-similar artificial networks can be 
classified as belonging to a unique universality class. 



B. Self-similar networks 

Let us now consider computer-generated models which 
satisfy Eq.([T|). For all these models, we find that Eq.® 
still holds. In strong contrast with non-self-similar net- 
works, the value of the scaling exponent i> now depends 
on the particular network analyzed and on the specific 
renormalization transformation. Moreover, each self- 
similar model is a fixed point of the renormalization flow 
since the statistical properties of the network are un- 
changed if iterated renormalization transformations are 
applied to the network. 
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FIG. 3: Study of renormalization flows on ER model with (fc) = 2. The box covering has been performed by using RB with 
r_g = 1 (a,b) and GCA with Ib = 1 (c,d). The figures display the susceptibility Xt (a,c) and the average clustering coefficient 
Ct (b,d) as a function of the relative network size xt. The insets display the scaling function of the variable [xt — x ')Ny" for 
Xt and Ct- Here the scaling exponent v = 2 and the susceptibility exponent "f — u in all cases. 



As a prototype of computer-generated self-similar net- 
work, we consider the Fractal Model (FM) introduced by 
Song et al. The FM is self-similar by design, as it is 
obtained by inverting the renormalization procedure. At 
each step, a node turns into a star, with a central hub and 
several nodes with degree one. Nodes of different stars 
can then be connected in two ways: with probability e 
one connects the hubs with each other, with probability 
1 — e a non-hub of a star is connected to a non-hub of 
the other. The resulting network is a tree with power 
law degree distribution, the exponent of which depends 
on the probability e. 

In the case of the FM network it is possible to de- 
rive the scaling exponent v, by inverting the construc- 
tion procedure of the graph. In this way one recovers 



graphs with identical structure at each renormalization 
step and one can predict how Kt, for instance, varies as 
the flow progresses. Since we are interested in renormal- 
izing the graph, our process is the time-reverse of the 
growth described in , and is characterized by the fol- 
lowing relations 

Nt-i = nNt, 
log s 

where n and s are time-independent constants determin- 
ing the value of the degree distribution exponent /? of 
the network. Here Nt and kt are the number of nodes 



6 




FIG. 4: Study of renormalization flows on the BA model with 2m — (k) — 6 (m indicates the number of connections introduced 
by each node during the construction of the BA model). The box covering has been performed by using RB with rs ~ 1 (a) 
and GCA with £b = 2 (b). The figures display the variables itt (a,b top) and rjt (a,b bottom) as a function of the relative 
network size xt. The insets display the scaling function of the variable {xt — x'')Nq^" for nt and 774. Here the scaling exponent 
i/ = 2 in both cases. 




FIG. 5: Study of renormalization fiows on the ER model with (fc) = 2. The box covering has been performed by using RB 
with rs = 1 (a) and GCA with =2 (b). The figures display the number of iterations t as a function of the relative network 
size Xt- The fixed point [i.e., Xt = x* = 0.059 (RB), 0.015 (GCA)] is reached in a number of renormalization steps growing 
logarithmically with the initial system size A^'o (see the insets on the right in each figure). In contrast, the number of stages 
needed to go out from the fixed-point (i.e., to reach a given xt < x*) grows as a power of A'o (see the insets on the left in each 
figure). 



and a characteristic degree of the network at step t of 
the renormalization; we choose the maximum degree Kt- 
The initial network has size A'o and shrinks due to box- 
covering transformations. In this case, for the variable 



Kt one obtains 
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FIG. 6: Statistical properties of the fixed point in the case of computer-generated networks. Renormalization has been performed 
by using GCA with £b = 2. Initial network sizes are: iVo = lO'* for the ER model, A'^o = 10® for the BA model, A'o = 10*^ for the 
Watts-Strogatz (WS) model (with ratio of rewired links p = 0.01), and A'o ~ 156251 for the fractal model (FM), with ratio of 
added connections p — 0.05. Dashed lines have slopes —1.5 in (a), —1.2 in (b) and —1 in (c). (d) The graphical representation 
of the fixed point has been obtained by starting from an ER model with Ao — 30000 and (k) = 2. 



where we used s = n^/^^ ^\ Ni/Nq — n * and Kq 



nI'^'^ ^\ derived from Eqs. 0. We see that the scaline; 
exponent v = 1 is obtained for any value of the expo- 
nent (3. From Eq.([5]) we actually get the fuU shape of 
the scaling function, that is a power law: our numerical 
calculations confirm this prediction (see Fig.[5)D). We re- 
mark that this holds only because one has used precisely 
the type of transformation that inverts the growth pro- 
cess of the fractal network. This amounts to applying 
the GCA with is = 3, as we did Fig. [Sb- 

If we consider instead the renormalization procedure 
defined by RB with tb = 1 (or by GCA with = 2), 
the centers of the boxes will be mostly low degree nodes, 
as discussed above. Hubs are thus included in boxes only 



as neighbors of low degree nodes and, as a consequence, 
the supernode corresponding to a box with a large hub 
inside will have a degree which is essentially the same 
as the degree of the hub before renormalization. It is 
therefore reasonable to assume that Kt ~ Kq, and we 
get 



Nt 



Nt 



N, 



i/(/3-i) 







|3-2 



Nt 



{Nr^t) 



(9) 



which is again a scaling function of the variable Ng^'^xt, 
with 1/ — (/3 — !)/(/? — 2), as we found numerically (see 
Fig.EJi). 

Qualitatively similar numerical results can be shown 
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FIG. 7: Study of renormalization flows on non-self-similar artiflcial graphs. Renormalization has been performed by using GCA 
with ^fl = 3. The figures display the variables Kt (a,b top) and rjt (a,b bottom) as a function of the relative network size xt, 
for the renormalization flow of an ER model with (k) = 2 (a) and a BA model with 2m = (k) = 6 (b). The insets display the 
scaling function of the variable {xt — x'')Nq^'^ for Kt and r]t- Here the scaling exponent is z/ = 2 in both cases. 




FIG. 8: Study of renormalization flows on self-similar artificial graphs. The box covering has been performed by using KB 
with rs = GCA with £b = 3 (b). The graph is an FM network with e = 0.5, where e is the probability for hub-hub 

attraction The figures display Kt (a, b, top), and xt (a, b, bottom) as a function of the relative network size xt. The scaling 
function of the variable {xt — x*)NIJ'^ for nt and xt is displayed in the insets. We find that the two box covering techniques 
yield different exponent values: v = 2.2 (RB) and v = 1 (GCA). The dashed lines indicate the predicted behavior of the scaling 
function. In (a) the exponent of the power law decay for the scaling function is —1, independently of the exponent /3 of the 
degree distribution of the initial graph; in (b) instead the scaling function decays with an exponent — (/3 — 2)/(/3 — 1) = —0.45. 
We still find ^ = v for both transformations. 



also for other self-similar models of networks: un- 
perturbed Watts- Strogatz (WS) model [jl (i.e., one- 
dimensional lattice), hierarchical model [23| and the 
Apollonian (AP) network model [l^l (see Table |T|. 



C. Effect of small perturbations on self-similar 
networks 

Self-similar objects correspond by definition to fixed 
points of the transformation. To study the nature of 
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Type 


Network 


K 


V 


X* 


Non-self-similar 


ER (fc> = 2 


Ib^2 
rB^2 


2.0(1) 
2.0(1) 
2.0(1) 
2.0(1) 


0.059(1) 
0.15(1) 






BA m — 3 


r_B = 1 
= 2 

rB=2 
=3 


2.0(1) 
2.0(1) 
2.0(1) 
2.0(1) 


0.098(2) 
0.245(5) 






Self-similar 


WS (fc> = 4 


rs = 1 
= 2 
^B =3 


1.0(1) 
1.0(1) 
1.0(1) 







FM e = 0.5 


rs = 1 
^B = 2 
^B =3 


2.2(1) 
2.2(1) 
1.0(1) 







AP 


^B = 2 
^B =3 


4.8(2) 
1.0(1) 






Perturbed self-similar 


WS (k) = 4 


rs = 1 
^B =3 


2.0(1) 
2.0(1) 


0.004(2) 



FM e = 0.5 


TB = 1 
^B =3 


2.1(1) 
2.0(1) 


0.118(2) 



AP 


rs = 1 
^B = 2 
^B =3 


2.0(1) 
2.0(1) 
2.0(1) 


0.045(2) 
0.05(1) 




TABLE I: We list the values of the scaling exponent v and 
of the fixed point threshold x* (fourth and fifth column, re- 
spectively) for all networks we consider in our numerical anal- 
ysis. Computer-generated networks (specified in the second 
column) are divided in non-self-similar, self-similar and per- 
turbed self-similar (first column). The perturbation is made 
by rewiring a fraction p = 0.01 of all links in the WS model 
and by adding a fraction p = 0.05 or p = 0.01 of all connec- 
tions in the FM or AP networks, respectively. The third col- 
umn specifies the type of transformation used to analyze the 
renormalization flow. We associate to each numerical value 
of V and x* its error. 



these fixed points, we have repeated the analysis of the 
renormalization flows for the self-similar networks con- 
sidered, but perturbed by a small amount of random- 
ness, through the addition or rewiring of a small frac- 
tion p of links. The results are shown in Fig. [5] for 
WS small- world networks, which are simply linear chains 
(trivially self-similar) perturbed by a certain amount of 
rewiring Q, and FM networks with randomly added 
links. In both cases we recover the behavior observed 
for non-self-similar graphs, with a scaling exponent v = 2 
(this holds for all values of rs or £b investigated, see also 
[201) ■ This clearly implies that the original self-similar 
fixed points are unstable with respect to disorder in the 
connections, and highlights once again the robustness of 
the exponent value v = 2. Furthermore, the statisti- 
cal properties of the fixed point reached at x*, when it 
exists (i.e., for = 1 or = 2) are again the same 
as those obtained starting from non-self-similar networks 
(see Fig. |6]). For these particular renormalization flows 
(for rB = lor^B = 2), the picture obtained is therefore 
a global flow towards the structure depicted in Fig. [SI 
with isolated unstable fixed points given by the artificial 



self-similar graphs. 



IV. REAL NETWORKS 

For real-world networks, a finite-size scaling analysis is 
not available because of the uniqueness of each sample. 
On the other hand, it is still possible to apply repeat- 
edly the renormalization transformation and to study the 
evolution of the network properties (a similar numerical 
study has been performed also in [29|]). In Fig. [TUl we 
measure some basic statistical properties of two real net- 
works along the renormalization flow. We consider the 
Actor Network (8i|ja graph constructed from the Internet 
Movie Database [S^ where nodes are connected if the 
corresponding actors were cast together in at least one 
movie, and the link graph of the Web pages of the domain 
of the University of Notre Dame (Indiana, USA) [26j . 
Both networks have been claimed to be self-similar, since 
Eq.([T]) holds for both of them On the one hand, 
the degree distributions P{k) are only slightly affected 
by the renormalization transformation, and retain their 
main characteristics even after several stages of renormal- 
ization (in particular for the Web graph, see Fig. [TOb). 
This first result points towards an effective self-similarity 
of P{k) under the action of the renormalization flow. 
The degree distribution by itself is however not enough 
to characterize complex networks, since many different 
topologies can correspond to the same P{k). Important 
information is in particular encoded in the clustering co- 
efficient spectrum C{k), defined as the average clustering 
coefficient of nodes of degree fc, and in the average de- 
gree of the neighbors of nodes of degree fc, fc„„(fc), which 
is a measure of the degree correlations between nearest 
neighbors in the graph. In this context. Fig. [TUb . d, e, 
and f show that even a single renormalization transfor- 
mation induces large changes in these quantities. In this 
respect, the apparent self-similarity exhibited by the de- 
gree distribution does not extend to higher order corre- 
lation patterns. 

Interestingly, in the case rs = 1 oi £b = 2, the it- 
eration of the renormalization transformation leads all 
real networks investigated [either self-similar or not, as 
defined by Eq. ([1])] towards the same kind of structure (il- 
lustrated in Fig. [6]) which is reached by non-self-similar 
artificial networks (see Fig. [TT|) . Note that, for real net- 
works, no change in the initial size can be performed, so 
we simply show the structure obtained after a few steps 
of renormalization, which remains stable for many steps 
due to the pecufiarity of the case rB = I or £b — 2, as 
explained above. 

All these results allow us to discuss the exact self- 
similarity of real- world networks: as we have seen in 
the case of computer-generated self-similar networks, 
all fixed-points correspond to strongly regular topolo- 
gies and minimal perturbations are enough to break the 
picture of self-similarity. Since randomness is an un- 
avoidable element in real complex networks, exact self- 
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FIG. 9: Effect of a small random perturbation on renormalization flows. The box covering has been performed by using RB 
with rs — 1. (a) WS network with (k) = 4 and a fraction p — 0.01 of randomly rewired links, (b) FM network with e = 0.5 
and a fraction p — 0.05 of added links. The flgures display Kt (a, b, top), and xt (a, b, bottom) as a function of the relative 
network size xt- Comparing with Fig.|8^, we see that the transformation yields a crossing of the Kj-curves for the FM networks. 
The crossing appears also for the WS networks, while in the unperturbed case the Kt-curves do not cross since these networks 
correspond to linear chains (trivially self-similar) . The exponents are now very different from the unperturbed case: we recover 
1/ = 2. The relation 7 = seems to hold here as well. 



similarity should not be observed in them. The random- 
ness of their topology is ampHfied when renormalization 
is iterated. Real-world networks which are self-similar 
according to Eq.(IT]) could however a priori be arbitrar- 
ily close to an actual fixed point of the renormalization. 
This actually raises the important issue of defining and 
measuring a distance in the space of networks. However 
Fig. [TO] shows that few renormalization steps (often a sin- 
gle one) are enough to substantially modify the network 
structure in our examples, so that it seems to us hard to 
sustain that they are close to a fixed point. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have presented a detailed analysis of 
the method, introduced in [2^, to study renormalization 
flows of complex networks. The method is applied using 
the two most popular techniques for network renormal- 
ization: greedy coloring [l^ and random burning [l^ . 
Independently of the algorithm, we have shown that a 
simple scaling rule [i.e., Eq.(I5])] is able to describe the 
renormalization flow of any computer-generated network. 
A single scaling exponent 1/ is needed in order to classify 
networks in universality classes: all non-self-similar net- 
works belong to the same universality class characterized 



hy ly = 2; self-similar networks, on the other hand, belong 
to other universality classes, generally identified by val- 
ues of the scaling exponent v different from 2. Self-similar 
networks represent by definition fixed points of the renor- 
malization transformation, since the statistical properties 
of these networks remain invariant after renormalization. 
They are actually unstable fixed points of the renormal- 
ization transformation: minimal random perturbations 
are indeed able to lead the fiow out of these fixed points. 
The numerical study presented here confirms and extends 
the validity of the results already anticipated in [2^ . 

In addition, we have performed an analysis of the ef- 
fect of the iterated renormalization transformation on 
real networks. Unfortunately, the same technique, in- 
troduced in [lOl, cannot be directly applied to real net- 
works. In this case, in fact, a finite-size scaling anal- 
ysis cannot be performed since any real network has a 
fixed size. Nevertheless, the usual simple measures of the 
network structure, taken after a few iterations of renor- 
malization, reveal that the transformation modifies the 
topological properties of the network. Furthermore, the 
repeated renormalization of real networks produces flows 
converging to the same fixed point structure, when it ex- 
ists, as the one found in the case of computer-generated 
non-self-similar networks. 
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FIG. 10: Statistical properties of real networks after t steps of renormalization. We consider two examples: a network of 
392, 340 actors, where nodes are connected if the corresponding actors were cast together in at least one movie (a, c and e); 
the link graph of the WWW, consisting of 325, 729 Web pages of the domain of the University of Notre Dame (Indiana, USA) 
and of their mutual hyperlinks '2^ (b, d and f). The box covering was performed with the GCA [Ib ~ 2), but the results 
hold as well for other transformations. The clustering spectrum and the degree correlation pattern change drastically already 
after a single transformation. In particular, the actor network displays assortativity, but after two transformations it becomes 
disassortative. The solid line in (b) has slope —2.1. 
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FIG. 11: Statistical properties of the 'fixed point' in the case of real networks. Renormalization has been performed by using 
GCA with £b ~ 2. The properties of the networks are measured after a certain number (less than ten) of renormalization 
steps. Dashed lines have the same slopes as those appearing in Fig. [S] The real networks considered in this figure are: the 
actor network the scientific collaboration network [s^, the network of Web pages of the domain of the University of 
Notre Dame [2a | and the protein-protein interaction network of the yeast Saccharomyces cerevisae [sj]. (d) The graphical 
representation of the fixed point has been obtained by starting from the protein-protein interaction network of the yeast. 
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